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Abstract 

We present a hyperbolic formulation of the evolution equations de- 
scribing non-radial perturbations of slowly rotating relativistic stars in 
the Regge- Wheeler gauge. We demonstrate the stability preperties of the 
new evolution set of equations and compute the polar w-modes for slowly 
rotating stars. 

1. Introduction 

The theory of non-radial perturbations of relativistic stars has been a 
field of intensive study for almost four decades, beginning with the pio- 
neering paper of Thorne & Campolattaro^ in the mid-60s followed by a 
sequence of papers by Thorne and collaborators. The authors focused 
mainly on perturbations of non-rotating relativistic stars. The study 
lasted for over three decades and today we are quite confident that we 
have a good understanding of the oscillation spectrum of non-rotating rel- 
ativistic starJ^. The study of rotating relativistic stars is a considerably 
more elaborate project. The construction of a rotating stellar background 
is by itself a very difficult task and during the last two decades a number of 
numerical codes were written for this purposJ^. The study of the response 
of a rotating star to small perturbations is also a very difficult analytic 
and numerical task. Fast rotating stars are not spherically symmetric ob- 
jects and the decomposition of the perturbations in spherical harmonics, 
which can in principle simplify the problem, is usually extremely difficult 
if not impossible. A way to study perturbations of rotating stars is to 
expand the perturbation equations in powers of the rotational parameter 
e = Q/Q,K, where Q is the angular velocity of the star and Ok is the 
angular velocity at the mass shedding limit. In this way, and up to 2nd 
order in e, the decomposition into spherical harmonics of the perturbation 
equations is feasible and the perturbation functions will be characterized 
by only one temporal and one spatial coordinate i.e. {t, r) . This "slow 
rotation approximation" is a very good tool in studying typical compact 
stars. For example, for the fastest known rotating pulsar with period of 
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1.56ms, we get e ~ 0.3. Yet compact stars during their birth could be 
rotating even faster and become dynamically or secularly unstable. For 
such cases and for the exact estimation of the critical angular velocity for 
the onset of these instabilities one might consider to do better than the 
"slow rotation" approximation and abolish the use of spherical harmonics. 
Then, due to the axisymmety of the background star, the perturbation 
functions will be characterized by two spatial coordinates (r, 9) and the 
temporal one t. Finally, with the advance of computational power it is 
now possible to evolve the full set of non-linear Einstein's equations de- 
scribing the rotating b ody an d from there to extract the characteristic 
frequencies of oscillatior l I I I Still in both cases i.e. the "fast rotation" 
and the "non-linear evolution" the extraction of the various features of 
the modes is quite cumbersome if not impossible without the proper un- 
derstanding of these features in the "slow rotation" approximation. This 
article serves this purpose, it presents a system of perturbation equations 
in the slow rotation approximation, which are stable and can be used to 
study the fundamental features of the rotating star perturbations. 

The full set of perturbation equations for rotating relativistic stars 
is quite complicated. During the last decade we got the first results for 
oscillation modes of rotating relativistic star j^ l ^ l ^ l ^^ l ^^ l ^'^ l '^^ l but there 
is still a long way to go before we have a full understanding of the be- 
haviour of the oscillation spectrum. In most studies mentioned earlier the 
slow-rotation approximation was used to tackle the problem. There are 
typically two methods to study the spectrum of the oscillation frequen- 
cies. The first method is by Fourier transforming the equations, where one 
drops the time dependence of the equations and the system can be solved 
as a boundary value problem, given that the correct boundary conditions 
have been provided. The second approach is the evolution of the sys- 
tem of linear (or non-linear) perturbation equations and then by Fourier 
transforming the perturbation functions one can get the characteristic fre- 
quencies which happen to be excited for a given set of initial data. Both 
techniques have their own advantages and disadvantages, for example in 
the first approach the correct definition of the boundary conditions is vital 
for the correct estimation of the frequencies, but this might be difficult for 
fast rotating compact stars. On the other hand via this method one can 
get the complete spectrum of oscillation modes. The time evolution tech- 
nique does not depend critically on accurate definition of the boundary 
conditions but the excitation of specific modes depends on the initial data 
provided. For a generic initial perturbation one should get all modes of 
the coupled spectrum, but they may be difficult to be found if excited at 
small amplitude or due to the finite evolution time. This means that one 
can reveal via this method only specific parts of the oscillation spectrum. 

Chandrasekhar & Ferrar^^ were the first to study axisymmetric per- 
turbations in the slow rotation approximation and apart from deriving the 
system of equations they have also showed how rotation induces coupling 
of the polar and axial modes, which actually decouple for non-rotating 
stars r Soon afterwards, Kojimsl^ presented the first complete derivation 



*PoIar or spheroidal or even parity modes are characterized by a sign change under parity 
transformation according to {—!)', while the axial or toroidal or odd parity ones change as 
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of the coupled polar and axial perturbation equations. But, in the specific 
form that they have been given, the equations where not extremely useful 
for time-evolutions. Kojima used the well-known Regge- Wheeler gauge 
for the spacetime perturbations. Allen et alM^ managed to derive from 
Kojima's set of equations a set of two coupled wave equations which to- 
gether with the appropriate constraints could descr ibe the perturbations 
of a non-rotating background. Finally, Ruoff et al.^l derived an alter- 
native set of evolution equations in a gauge suggest ed t hirty years ago 
by Battiston et aJ^. In their derivation Ruoff et alM^ used the ADM 
form of linearized Einstein's equation^^ and showed that the use of the 
Regge- Wheeler gauge leads to unnecessary terms with second order spa- 
tial derivatives of the lapse function a and the 3-spatial metric 7ij. The 
Battiston et.al. (BCL) gauge, appears to be more suitable for the lin- 
earized ADM formulation of Einstein's equations but it seems to suffer 
from numerical instabilities which have not yet been resolved. It is not 
actually clear whether the numerical instabilities observed are due to the 
known instability of the ADM formulation of Einstein's equatio ns o r it is 
due to a combination of the formulation and the choice of gang 

This is the reason that we derived a system of hyperbolic equations 
describing the linear perturbations of slowly rotating stars using the stan- 
dard form of Einstein's equations (not the ADM formulation) in the com- 
monly used Regge- Wheeler gauge. These equations reduce, in the non- 
rotating case, directly to the equations of Allen et al!^ (polar pertur- 
bations) and in the slow-rotating case to the axial ones that have been 
used by Anderssort — ' and Ruoff & Kokkotas — '-^ in the study of the 
relativistic r-mode spectrum. 

As a test of stability of the present system, we have derived a reduced 
version of the equations describing only the space time perturbations in 
the so called inverse Cowling Approximation (ICA]p31. The fluid part is 
actually unaffected by the choice of gauge. We show that in the present 
form the system of evolution equations is stable. Additionally, we have 
calculated, for the first time, the polar w-modes for rotating stars. 

In the next parts of this section we briefly describe the equations for 
the construction of the background stellar model and the derivation of 
the perturbation equations. In the second section, we show the form of 
the evolution equations for both fluid and spacetime perturbations. The 
derivation of the evolution equations for the fluid perturbations, in a for- 
mat compatible with the equations describing the spacetime perturbations 
is given in Appendix I. For clarity, we have chosen to give the lengthy 
background coefficients of the evolution equations in Appendix II at the 
end. In the third section, we describe a simplified form of the evolution 
equations for the spacetime, i.e. we "freeze" the fluid perturbations and 
keep only the part of the system that deals with the spacetime equations 
(ICA). We evolve this system of equations as a test of stability and com- 
pute the w-modeJ^ for some typical rotating stellar models. Finally, in 
the last section we discuss the analytic and numerical results that have 
been presented. 
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1.1. The background model 

Here we will briefly present the basic equations describing the station- 
ary configuration. We consider a uniformly rotating relativistic star with 
angular velocity fl, which is assumed to be small compared to Qk, the 
angular velocity of the mass shedding limit. We will keep terms of first or- 
der in Q only. The parameter that we will use to expand all the rotational 
quantities will be e = H/^Ik- Higher order terms in the construction of 
the background models as well as in the perturbation equations will be 
omitted. 

Under the above assumptions the background metric of a slowly ro- 
tating spherically symmetric star has the form 

= -e"dt^ + e^dr'^ + r\de^ + sin Bdcj)^) - Iwr^ sin^ edtd(i>, (1) 

where v and A are functions of the radial coordinate r. The function 
a; = aj(r) describes the dragging of the inertial frames and is a first-order 
rotational correction. 

We will also assume a perfect fluid described by the energy-momentum 
tensor: 

Ty.y = (p 4- p)u^Uy -I- pg^y (2) 
with a 4-velocity vector of the form 

= (e-"'/^ 0, 0, J^e-/") (3) 

where p and p are the pressure and the energy density. The static spherical 
stellar conflguration will be constructed using the Tolman-Oppenheimer- 
Volkof (TOV) equations'^^ . The solution of the TOV equations provides 
the metric functions (potentials) v = i/(r) and A = A(r) and the remaining 
function = w(r) (dragging) will be calculated in the stellar interior via 
the equatiorP^ 

lJ' — [k(p + p)re^ -f 4r^^] u — 4k{p + p)e'^uj — (4) 

where k = 4it. Outside the star, lo gets the simple form ui — 2J/r^, where 
J is the angular momentum of the star. An additional function 

zu = O, — uj . (5) 

will be used in various equations for simplicity reasons. Note, that in the 
present study we will use a barotropic form of the equation of state i.e. 
p = p{p) , which serves well for the problems that we consider here. 

1.2. The perturbed configuration 

For the derivation of the perturbation equations we assume a a metric 
perturbation fe^^ which is small compared to the background metric g'^J 
and the perturbed spacetime will be described by a metric of the form 

Sm-' = glS + ht,y . (6) 
Still we will use for lowering and raising the indices the background metric 
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For the fluid perturbations we assume a small velocity perturbation 
5u'^ and small density {5p) and pressure {5p) variations. With these as- 
sumptions the perturbed Einstein equations 

<5Gp. = 2kST^^ (7) 

together with the equation of motion for the fluid perturbations 

<5(T''^^)=0 (8) 

provide a complete set of equations for the study of small fluid and space- 
time perturbations in a spherically symmetric, slowly and uniformly ro- 
tating background. 

The components of the perturbation tensor /i^j, will be functions of 
all four coordinates {t, r, 9, 0) but due to the spherical symmetry one can 
expand its components in tensor spherical harmonics. Then by making 
an appropriate choice of gauge, here we choose the Regge- Wheeler one, 
ftpi/ gets the form: 



oo L 



( e"Ho,im Hi^im -/lo.imsin '^684, ho^imsinede 



i = 2 m = — l 



r'^Ki^ 



\ * * * sin^ dKim 

(9) 

where Ho{t,r), Hi{t,r), H2{t,r) and K{t,r) are the functions describing 
the polar perturbations, ho{t,r), h\{t,r) describe the axial ones and the 
star symbol (*) stands for the spherically symmetric comp onen ts of h^v 
The contravariant form of the perturbed fluid velocity i 



Su = 



; = 2 m. = -l 



"'^Ho^im + ^e '^"^^ ho^im sin ede^ Yi^ 

1 — 2 m— — l 
,./2-A 

r r ^ 

= 77 7 

k(p + p) 

^ fc J+ E E (^""^'^ - ^^^) (12) 

i,/2 °° ' 

^""^ = 77 r-r^^y y (V^im90 + (7,„ sin eae)y,™ (13) 



1 = 2 m=-l 



1=2 m = — l 



where again R{t,r), V{t,r) are polar perturbation functions and U{t,r) 
describes the axial fluid perturbation. Finally, for the density and pressure 
perturbations we will use the conventions 
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For the barotropic case that we study, the sound speed relates the varia- 
tions of pressure and density i.e. 

Concluding, our task will be to derive a system of evolution equations 
which will involve the five spacetime perturbation functions (Hq, Hi,H2, K)i,-a, 
{ho, hi)im, the three velocity perturbations {R,V,U)im and the density 
perturbation pi^im- There are various ways of writing the system, for ex- 
ample one can derive a system of evolution equations for every single one 
of the above functions. Another way will be to evolve a certain number 
of the functions and derive the rest from the constraint equations. Here 
we choose the first approach since at least in the non-rotating case it has 
been proven to be numerically more 

stabl«Pl. 

2. Perturbation equations 

The generic form of the perturbed Einstein equations in the back- 
ground of a slowly rotati ng star i.e. the analytic form of equation Q is 
actually given by Kojima-^. We derive in Appendix I the equations for 
the fiuid perturbations i.e. from equation ||5J which were not derived by 
Kojima. Then we perform the appropriate transformations in order to 
derive a system of first order evolution equations for the spacetime and 
fluid perturbation functions. 

The difficulties in the derivation of the evolution equations and in 
their numerical evolution comes from the rotation effects which dramati- 
cally affect the structure of the system and as a consequence the oscillation 
spectra. Rotation affects the spectrum in various ways. First, it erases 
the azimuthal degeneracy with respect to m, second it adds rotational 
corrections to the equations describing perturbations of non-rotating con- 
figuration and third it couples the perturbation equations of a specific 
harmonic index I with the / ± 1 perturbations, i.e. it couples the polar to 
the axial perturbations and vice-versa. In this way the quadrupole {I = 2) 
polar (spheroidal or even parity) perturbations couple with the octupole 
{I = 3) and dipole (/ = 1) axial (toroidal or odd parity) perturbations! 
Therefore the system of evolution equations, keeping rotational terms up 
to 1st order in e and for a specific harmonic index I will be of the following 
form: 



Pirn = {&'°'+&'^')PL+{p'p°'+P'p'')Pl 



Im 



— (Q^A +'^A ]^lm+\PA +Pa ] ^Im 

+ 7« {cf,Ct,Ct) PL + 5^p^ {ct,Ct,Ct) Pirn (17) 



where the dot corresponds to the time derivative and the prime to the 
spatial derivative with respect to the radial coordinate. The arrays Pim = 



^This is a true if we keep rotational terms up to first order in e, the inclusion of higher order 
terms makes the coupUng more involved i.e. it adds coupling terms of i ± 2 
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{Kim, Ho,lm, Hl^lm, H2,lm, Rlm,Vlm) and Aim ~ {hi ^Im , ho,lm , Ulm) Stand 

for the functions which describe polar and axial perturbations correspond- 
ingly. The coefficient matrices of the form Op ' , d^' , /3p ' and Z?^' include 
only background terms without rotational corrections, while the matrices 

id) ^(1) 

A,W A,W xW _j iW 



of the form dp', d p'p^ and P'p'' are the ones which include the rota- 
tional corrections. Finally, the matrices 7p', 7^^', (Jp' and iJp' include 



background rotational terms and especially the operators Cf , and £, 
which are responsible for the coupling with the I ± 1 perturbations. The 
analytic form of these three operators i 

L^Alm ~ {I ~ l)QlmAl-l m — {I + 2)Q; + l rnAl + l m (18) 

Aim ~ —{I + i)QlmAl-i m + lQl + 1 m^! + l m (19) 

/Is Aim = {l + l){l-l)QlmAl-lm+l{l + 2)Ql + 1 + 1 m (20) 

with 



{l-m){l + m) 



{21 - 1){21 + 1) 

We should stress that because of the coupling between the Z-polar per- 
turbations with the I ± 1-axial and vice versa the above sets of equations 
lllbll - 117|l have to be considered as infinite system with / running from \m\ 
to infinity, when m is considered to be fixed. There is no coupling to equa- 
tions with / < |m|, since from its definition (I21II it follows that the relevant 
couphng coefficient Qmm is zero. Furthermore, the equations form two 
independent sets, each belonging to a different parity. In the non-rotating 
case one usually distinguishes between polar and axial perturbations. In 
the rotating case, however, the polar and axial equations are coupled, and 
therefore the distinction between polar and axial modes can no longer be 
upheld. However, the equations do not mix the overall parity, which can 
be seen as follows. A polar equation with even / — \m\ has even parity. It 
is coupled to an axial equation with 1 + 1, whose parity is also even. This 
axial equation in turn couples to a polar equation of order 1 + 2, which, 
again, has even parity. As this continues in the same manner, we can 
see that for even m, the complete coupled system with a leading polar 
equation has even parity. Conversely, the other system starting with a 
leading axial equation has odd parity. For odd m we obtain the reversed 

situation. This implies that for any given m, it makes sense to distinguish 

6 

the modes according to their overall parity. Lockitch and Friedman in- 
troduced the notion of polar or axial led modes, depending on whether 
the leading equations with I = \m\ are axial or polar. For even m the 
polar led modes have even parity and the axial modes odd parity, and 
vice versa for odd m. The axisymmetric case m = is somewhat special 
in that only the even parity modes start with I — whereas the odd parity 
modes have to start with / = 1, for there are no axial I = equations. 

In the derivation of the wave equations describing the spheroidal (po- 
lar) oscillations of non-rotating stars Allen et all^^had derived a set of two 
wave equations for the following combination of the metric perturbations 

Flm=rKlm, Sim = {Ho,lm ~ K^n) (22) 
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while the enthalpy 

Him = (23) 
P + P 

was calculated for every time step by solving an elliptic equation, the 
Hamiltonian constraint. They have observed that the evolution was more 
stable when they were using evolution equations for all the perturbation 
functions F, S and H. This is the reason that in our derivation of the 
equations we choose to produce an system of evolution equations for all 
the function describing both the fluid and the spacetime perturbations. A 
similar set of evolution equations for non-rotating stars derived recently 
by Nagar et.al.^^ while for gauge invariant formulations one can refer to 
a series of articles by Gundlach and Martin-Garchia' '^ I '^^ i '^ i 

Using a similar approach to the one used by Allen et al!^ we derived 
a set of wave equations for the rotational case. The form of the evolution 
equations for the spacetime perturbations is: 

F-e^-^F" = ao{F,S,H,F')+imai{R,V,HuV',H[,F,S) 

+ Cta2{ho,hi,U,h'o,h[,U',hi) (24) 

S-e^-^S" = boiS,F,S')+imbi{R,V,Hi,R',V',H,F,S,V) 

+ Cfb2[ho,hi,U,h'Q,h'i,U',ho,hi,h'Q,h'i) (25) 

Hi = co{S,F,S')+imci{R,V,Hi,F,F') 

+ Cfc2{ho,U,h'o,ho,hi) (26) 

'h\~e''~^h'l — do{hi,h'i) + imdi (^ho,hi,U,h'(),h'i,U' ,hij 

+ Ctd2{R,V,Hi,V',F) (27) 

ho = eo{hi,h'i) + imei (^ho,U,h'o,hi) 

+ Cfd2{V,F) (28) 

Where we have used the following rule in naming the various functional 
forms which are involved in the above set of equations. The ones of the 
form ao, foo, co, do and eo do not include rotational corrections i.e. the 
coefficients of the perturbations functions are of zeroth order background 
terms. Those of the form ai, b\, ci, d\ and ei include rotational cor- 
rections and related to the perturbation functions which have the same 
parity as the evolution function. Finally, those of the form 02, 62, C2, d2, 
62 introduce rotational corrections and couple with the perturbation func- 
tions of the other parity i.e. with spherical harmonic indices Z ± 1. This 
last set of coefficients includes the operators , and CL^ which are 
the ones that introduce the coupling with the Z±l perturbation functions. 



2.1. The first order system of equations 

It is evident that the above system of equation is not in a form suitable 
for time evolutions since a significant number of equations involve time 
derivatives or even mixed derivatives on the right hand side. These terms 
are due to the rotational corrections of the equations and they were not 
present in the non-rotational case. We will abandon this writing (2nd 
order in space and time wherever possible) and instead we will derive a 
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first order in time hyperbolic system of equations. Still the above sys- 
tem of equations is useful in its present form for time evolutions of the 
perturbations of non-rotating relativistic stars. 
We proceed by defining the new functions 

Fi =F, Si = S, Ti = fti 

F2 = F', S2 = S', T2 = h[ (29) 

This additional set of functions combined with appropriate substitutions 
provides the new set of evolution equations. The equations describing the 
evolution of the spacetime perturbations will be: 

Fi = e'"^F2+aoiF + ao2S + ao3H + ao4F2 

+ aio-R + aiiV + ai2-ffi + aiaV + anH[ + ais-Fi + aisft 

+ a2o/io + 021/11 + a22?7 + a23/io + a24/i'i + a25?7' + 026^1 (30) 



Si = e-''^S2+boiS + bo2F + bo3S2 

+ bioR + buV + 612H1 + biaR' + biiV' + bizH[ + bwFi + b^rSi 

+ 620/10 + &21/II + 622?/ + 623/I0 + 624/I1 + &26?7' + 626Ti 

+ b27Ti + 628T2 (31) 

Hi = cooF + coiS + C02S2 

+ cio-ffi + ciiV + C12R + C13F1 + C14F1' 

+ C20/IO + 021^ + C22/lo + C23Ti + C24/I1 + C26/ll (32) 

Ti = e'~^T2 + doi/ii + do2T2 

+ dio/io + dii/ii + di2U + di^h'o + di4T2 + disfj' + cJisTi 

+ d2oHi + d2iV + d22-R + d23V' + d24-Fi (33) 

ho = emhi + eoi/ii + eiofto + euU + 612/10 + eisTi 

+ eaoV + e2iFi (34) 

The derivation of the fluid perturbation equations is presented in Ap- 
pendix A. The equations derived there can be combined appropriately 
with the above system equations to provide a system of four evolution 
equations for the perturbation of the enthalpy function i?;^ and the three 
velocity perturbation functions Rim, Vim and Uim- They actually get the 
following form: 

H = koiR + fco2 V + koaFi + hnSi + hK>H[ + koeR' + koiHi 
+ kioH + kiiS + ki2F + ki3F2 + kiiFi 

+ k2ihi + k22T2 (35) 



R = fooF' + foiS' + fmH' + foaF + foiS 

+ fioHi +fiiV + fi2R + fisFi + hiF[ + fisV 
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+ /20 h'o + /21 ho + f22U + /23 U' + /24/ii + /25 h[ + /aeTi (36) 

V = qooH + qoiF + qo2S 

+ qioHi + quH'i + qi2V + qnR + qi4.Fi + qisSi + qmR' 

+ q2oU + q2iho (37) 

U = Soo/ll + Soi/l'i 

+ sio?7 + sii/io + S12/10 + S13T1 

+ S20V + S2lR+S22Hi+S23H[+S2iR' + S25F1+S26S1 (38) 

Here the various coefficients involved have the meaning described earher. 
In other words the ones of the form aoi, &oi, coi, doi, eoi, koi, foi, qoi and soi 
do not include rotational corrections. Those of the form an, bn, cu, du, 
Cii, kii, fii, qii and su include rotational corrections and related to the 
perturbation functions which have the same parity as the function that is 
evolved in every equation. Finally, those of the form a2i, b2i, C2i, d2i, e2i, 
k2i, f2i, q2i and S2i include also rotational corrections and correspond to 
the other parity and for spherical harmonic indices / ± 1. This last set of 
coefficients includes the operators jCf , and which are responsible 
for the coupling with the / ± 1 perturbations. All these coefficients are 
given explicitly in Appendix II. 

The above system of evolution equation is, as exp ecte d, considerably 
more complicated than the one derived by RuofT et alP^, using the BCL 
gauge in the ADM form of Einstein's equations. Still, it has considerable 
merits since parts of this system have been tested and show to be sta- 
ble. For example, in the non-rotational case it has been used by Allen 
et alM^. The part that corresponds to ax;ial perturbations (includi ng ro- 
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tations) have been use d in earlier papers by Ruoff and Kokkotas ' ' 
and Kokkotas et al.^^, w hile the evolution equations for the fluid have 
been tested by Ruoff et al.l^ in the Cowling approximation. Although, 
we don't use in this paper the above system of equations in its full form, 
we will test, in the next section, the part that deals with the evolution of 
spacetime perturbations. This part is shown to be numerically unstable in 
the BCL gauge, while in the present form it shows stability which allows 
us to calculate the w-modes for rotating relativistic stars. 

3. Inverse Cowling Approximation 

As we discussed earlier a flrst and crucial test of the evolution equa- 
tions derived here will be the evolution of the spacetime equations decou- 
pled from the fluid ones. This simplification is named Inverse Cowlirii 
Approximation (ICA) and it has been suggested by Andersson et al! 
as an easier way to understand t he spacetime perturba tioiis a nd their 
characteristic modes, the w-modej^. It has been showrJMl that the 
w-modes do not excite significant fluid motions and the omission of the 
fluid perturbations is a legitir nate approximation. 

As discussed by Ruoff et alf^the various families of fluid modes which 
are not degenerate in the non-rotating limit i.e. the f, p and g-modes are 
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not afTected in a crucial way by the coupling to the / ± 1 perturbations. 
This is not true for the inertial (or rotational) modes, for those the cou- 
pling to the various I ± 1 terms is vital for a proper understandin g of t he 
spectrum, see also a detailed discussion by Lockitch and FriedmarP^lsJ 

The w-modes are similar to the quasi-normal modes of the black holes 
where the rotation splits and shifts the spectra but (at least in slow- 
rotation limit) the coupling with the various Vs does not affect in a critical 
way the spectrum. Based on the above arguments we studied the w-mode 
spectrum of rotating relativistic stars. 

The simplified version of the evolution equations for the spacetime 
perturbations is: 

F = Fx F2=Fi (39) 
S = Si S2 = S'i (40) 

Fi = e^^^Fa + ^ {3/3 + a + 4M - Ar) F 

+ '^{a + I3~re-^) S +^{a- f3)F2 

+ i!^[(a-/3)Hi+rV^//(-rV^5i] 

+ ^ [(q + /3 - rA - re~^)t^ - r^e" V] Fi (41) 

Si = e''-^S2 + ^{3a + l3~rA)S 

+ ^(e^a^-2Mr-f-r/3)F-f^(a-/3)S2 

_ 2im^.-x^,jj^ _ 2im^ . + /J) + Ar - 2re-^l Si 
r rA 

+ ^^^[2r'e-^uj' + 2zu{a + P)+u;{a~f3~2R4)]Fi (42) 

Hi = ^e-'+^F + S + rS2 

- imuHi + ^ [(206^ + r)uj - r'J\ Fi + '^r^cjFl (43) 
rA J rA 



where 



A = ;(; -f 1) , a = M + kpr^ and (3 = - [M - kpr^) (44) 



Additionally, the following Hamiltonian constraint equation can be used 
to set up the initial data for the evolution: 

Fa - ^eVa + ^ (3/3 + 2A4 - rA) F 

+ — e^"" [4(0 + 13)- 4re"^ - rAl 5* - re'" S2 
2r 

+ im^A-. r^^^^ ^ ^ _ ^g-A^ _ 3^^^ _ 2r^e"^aj'l FI 
rA 

— imtjr e oi + — — e 2uj{ol — I3)e — r uj \ Hi 

2r -' 

— "^^^e'" ruoF'i + imijjre~" H'l = (45) 
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e 


n 










0.0 





10.8040 


10.8040 


10.8040 


10.8040 


0.2 


0.4995 


11.5294 


11.1413 


10.5070 


10.2247 


0.4 


0.9902 


12.5044 


11.5294 


10.2245 


9.7312 


0.6 


1.4853 


13.4466 


11.9802 


9.9763 


9.3629 


0.8 


1.9804 


14.2875 


12.4622 


9.7311 


8.9827 


1.0 


2.4755 


14.9495 


12.9565 


9.5126 


8.7620 


1.2 


2.9706 


15.5042 


13.4471 


9.3263 


8.5791 



Table 1: The frequencies of the I — 2 polar w-mode for various rotation rates 
for the polytropic model i?=8.86km and M = 1.267Mq. The frequencies of the 
modes and the rotational velocity are given in kHz. 



In our numerical study we have used polytropic equation of state of 
the form p — Kp^ and uniform density stellar models with typical masses 
and radii. The numerical code was stable for all values of the rotational 
parameter e. 

In Fig. I we present the evolution of the frequency of the first w-mode 
as a function of the rotational parameter e. The model that we have used 
has a polytropic equation of state with V = 2, K = lOOkm"^, radius R — 
8.86km and mass M = 1.27M0. We observe the splitting of the frequency 
for m = ±1 and m = ±2, while the frequency increases/decreases linearly 
with respect to the rotational parameter e. In the level of approximation 
that we have used the axisymmetric modes (m = 0) are not affected by 
rotation. 

In general the frequencies of rotating stars follow a pattern of the form 

o" = (To — KmQ, (46) 

where ctq is the frequency of the non-rotating star and «: = k(M, R, ...) is 
a function depending on the details of the specific stellar model. 

In Fig. 2 we show the frequency for two uniform density models. Both 
have mass M = 1.4Mq one with M/R = 0.20 and a more compact one 
with M/R — 0.25. One can easily observe that the more compact the 
model the higher is the w-mode frequency ctq of the non-rotating star as 
was shown by Kokkotas and Schuta— ^ 

The rotational corrections with m > for high rotation rates can lead 
to negative frequencies of the modes which signal the onset of an insta- 
bility. It is apparent from Fig. 2 that the more compact the model is the 
faster the frequency tends to zero i.e. to the instability. For the models 
discussed here the instability will never be excited for physically accept- 
able values of rotation i.e. for e < I. Axial w-modes for the polytropic 
model used here a s well as for uniform density stars have been calcu- 
lated alreadjE^El. Actually, the ultra-compact uniform density models, 
M/R > 0.4, can become secularly unstable as has been shown by Kokko- 
tas et alP2. 
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1.2 



Fig. 1. The 1 = 2 polar w-mode frequencies of a rotating polytropic star (R 
M = 1.267Mq) as functions of the rotation parameter e. 



8.86km and 



4. Discussion 

In this article we have studied perturbation of slowly rotating relativis- 
tic stars. In an attempt to find a stable system of evolution equations we 
have derived a quite long but seemingly stable hyperbolic system. In our 
calculations we have rederived the linearized form of Einstein's equati ons 
for a slowly rotating star. This form is identical to the one by KojimcPl 
and for this purpose we refer the reader to his article. For the evolu- 
tion of the fluid perturbations we have derived the full set of equations in 
Appendix I. 

Finally, we have tested a reduced version of the above system by cal- 
culating the w-mode frequencies for slowly rotating stars in the Inverse 
Cowling Approximation (ICA). 

Concluding we have to admit that although the various parts of the 
derived hyperbolic system of equations are well tested and are stable, the 
whole system of equations involving the coupling with higher and lower 
I's is yet to be studied. This by itself is a difficult task though the benefits 
suggest that it worth the effort. An interesting question that has to be 
answered is whether the instability observed in the ADM formalism of 
the equations in the BCL gauge is due to the choice of the gauge or the 
fundamental instability of the ADM formalism. 

From the full system of equations one can study all types of modes (f, p, 
w and inertial) including their rotational corrections and their couplings. 
This is a very useful information in understanding the behavior of the 
various modes as function of the rotation of the star. It can actually 
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Fig. 2. The I = 2 polar w-mode frequencies of two uniform density rotating stars {R = 10.34fcm 
and M = IAMq) as functions of the rotation parameter e. Both models have mass M = 1.4M0, 
the more comact model has radius R = 8.27km (M/R = 0.25) and the other R = 10.34km 
(M/R = 0.20). 



contribute in the gravitational asteroseismology '"^ * ^ ' ^ * by predicting the 
angular velocity of the star via the shifting of the various mode frequencies. 
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5. Appendix I: Equations for fluid perturbations 

In this Appendix we derive the perturbations of equations of motion 
for the fluid. The will be derived from the variation of the conservations 
laws of the energy-momentum tensor, i.e. from equation 10: 



5(r" 



= 



(47) 



The variation of the equations describing the conservation of energy and 
the radial component of the momentum will be written as: 



(48) 



where 7 = 4, 5. We use these values for / because t he a bove two equations 
can have the same form as equation (14) in KojimiP^ and can be consid- 
ered as a complementary set to these equations. Next we shall decompose 
the above equations to a specific mode with I, m. Multiplying equation 
Ij48l with Yi*^ and integrating over the solid angle, we get; 



A\2Yim + sm{e) + Qi^,^{l - 1)B 



I — Im 



-Q,+i^(Z + 2)bW,„ 



which has similar structure as equation (20) of KojimsPI. In 
tion this equation will be written as: 





(49) 
our conven- 



A«y,„ + \mC\2 sm{d) + CfB^^ = 



(50) 



The variations of the fluid equations for the {6) and {<j>) components 
have the form: 



+ 

+ 



^im' sin dYl^ + ifraXlm 

x['2 sm(e)H/,„ = 



edYr 



(51) 



im + cos(6») j sin ^ 9dYi^ 

Cim' sin OYlrn + XlraXlm 

- i[2^Me)Wim = Q (52) 

which have the structure of equations (15) and (16) of Kojima^, there- 
fore we will index them as J = 2. Their decomposition into spherical 
harmonics leads to equations of the form 



(53) 







im 


"(A- 






AJ)' 










r± 

'-2 


(J) 


-(A± 




+ ^3 «L 


= 






im 


"(A- 




+ 42 + 


'/im 










r± 

'-2 


'AJ) 


+ (A± 


'^12 


+ c 


±3(J) 

3 Plm 


= 



(54) 



where A — l{l + 1) and A* stands for A where l^l + lorl^l — 1 i.e. 
A+ — [l + 1){1 + 2) and A^ ~ {I — 1)^- T he a bove set of equations has the 
same structure as equations (24)-(25) in^l 
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The various functionals used in the above sets of equations have the 
following form: 

= -pi + 2pHo-U-'R' + e-\p-p)H[+'^V 



3(4) 



^1 



(4) 
Im 



2kr 



kr'^ 

(3ri^' - r\' + 4)R 



ip + p) (k+^H2 



%^ {rv'{^ + p)C7' - 2rv'p + (p - p)(rA' - 4)] i/i 



2r 



(55) 



{(p + p) [(n(rA'+rC7^!^'-4) -2ra>(A' + z/')] +8a;p}/ii 



2 (fip + vsp) K 



-(p + p)zu + 2pL0 



Ho — piO - piw 



(56) 



(57) 



4(5) _ 
■^Im — 



,(5) 



+ 



r^(5) 



+ 2e- "pi/i + 2pi/^ - -pK - u'pHo + p[ 
k r 

+ y (Pi + {^P - pv') H2-^ (p + p) i^o 

26^" (pfi + pzu) ill + ^ irwv' + rw' — 2zu] U 
kr ^ ' 

e-'{p + p)nh'o 

2 1 
{p - p){lo' - ojv') - -{p -\- p)w ho 



(58) 



(59) 



\q,R + 26"" (Op + rop) + ^ (rz^'ro + ruj' - 2w) V (60) 



(2) 
Zm 



3(2) 
(2) 

^(2) 

(2) 

a(2) 



fee"" [2wp - (p + p)n] /lo + (2w - a)f/ 
f/ + 2ke~'"pho + 2kpe~^h[ - ke~^ (pX' + pv' - ^p 

(2w - n)v 

fce"" A (^^{p + p)Q-pu?jho + ^AnU 



(61) 
(62) 
^hi (63) 

(64) 
(65) 



= ke (^^{p + p)Q. - pujj ho + ^iiU 



kr^e-^- 



--{(p + p) p - - 1) ^ - 1 (C77^z^' - A') 
+ w (^pA' + pi^'+ ^(p-p)^ 

+ ire"^ [4(0 + w) - 2rw' + (2a) + ft) ri/' - A'rO] 



(66) 
(67) 
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- ^kr'^e " \{p + p)vj - ^puj\ Ho + kr^e " {wpi + fipi) 

+ kr^e'" [3{p + p)Q. - 2puj] K 

+ fcr^e"^-" [(p + p)w + (p - p)a;] H[ . 



(68) 



Some further simplifications are still possible in the above two rela- 
tions. Therefore the terms related to the rotational corrections get the 
form: 



^3 O^lm 



and 



+ £-2 



(J) 



(.J) 



= ke^"" {[n{p - p) - 2zup] Cf + [(p + p)(A^ + l)n-2pa)]£^}/io 
+ [{n - 2vj)£.f + nC^] U (69) 



(J) 



2pu; - pKoj + -(3A - 2)(p + p) 



ho 



+ i(A-2)ny-|-i[4a; + (A-2)n][7 .(70) 

In section 2 the above set of equations is simplified still further and 
we will use instead of the functions Kim and i?o,im the new ones i.e. Fim 
and Sim- 
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6. Appendix II: Coefficients of the evolution equations 

In this appendix we provide the background coefficients for the evolu- 
tion equations l|29|l - Ij88|l . 

We have used the following substitutions for the pressure and density 

a = M + kpr'-^ = \r'^e''^v' and (3 = - (M - kpr^) = ^^^e^^A' (71) 

which simplify considerably some of the equations. 

Background coefficients for the evolution equation 13UII for Fi 

aoi = ^[3/3 + a + 4M-Ar] (72) 

ao2 = '^{a + P-re-'') (73) 

ao3 = ^e''(a + /3)(C7' - 1)) (74) 

ao4 = -j(a-/3) (75) 
Sim 2 

aio — — r vj (7d) 

^ , 16im „ / Sim 2 / 
flu = 4i'mre ui H — e cc7 (^p — re j r tz? (77) 

ai2 = — (a-/3)t^ (78) 
r 

Sim 2 .-rix 

ai3 = r ~ — aio (79) 

ai4 = imroje"'^ (80) 

ai5 . 

rA 

ai6 = — imr e^"u (82) 

4 2 
a20 = ^ [ai^ + P'^] -Ci + -7-^ ('ir^e~^Lu' — 2acj) £3 



"''"^ [(a + - rA - re-^)a; - r'e-^'aj'] (81) 



r 



Ar2 



+ ^e^^o;' [Sre^^ + 4(2q + 3/3) + r(A± - A - 1)] £^ 

+ ^ (q + /3) [2qC7' - Ar + 2(3q + 2/3)] tzj^J 

+ }^ [A*(^e-^-/3)-a]a;£± (83) 

"21 = ^;^(a-/3) (^2/:f +3^£±^ (84) 

a22 = 4re''ci7£f + — [(4/3-4re"^ -rA)ti7 + 4e~^r^a;'] £j (85) 

023 = -— [2r(l-A*)a;-16ti7(Q + /3)-(8re"^-rA + 4/3)ra;'] £^ 
- 2e'^{Lo + rLu')Ct + '^e'^LoCt (86) 

024 = ue''-^'' (^2Cf +i^Ct^ (87) 
r± 

0.25 = — Se r zu- 



A 
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-A I r± I 2 -A r± 
026 = re Lo + -^e loL^ 

-A 

+ — [8(a + f3)uj + r(8re"^ - r A + 4/?)a;' - r(A* - 2)lj] £^^=(89) 
Background coefScients for the evolution equation for 5*1 



&01 



6o3 = 
bio = 
+ 



6l3 = 

bi4 = 

615 = 

616 = 

+ 



+ 



-- (3q + /3 - rA) 



feo2 = — g- (e^Q^ - 2Mr + r/3) 



4i"i 21-. , o N/o a , o -An^2 , 8im 2^-a / 

— e (tj + 2ia7)(3q — p + 2r-e jCs H r—re vj 

rA A 



4im 
rA 
Sim 



e^" [2c:j(2a + re~^ - /3) + (3a - /3 + 2re"^)Lj] 



6ii = ~TT^^'^ [(20^6"^ + 6ra — Ar^) + (q + re ro'] 



r^A 
4im 

r 
4im 



612 = ^e''Cl{2vj + Lo)[{a-l3~ 2M) (2re-^ - /?) - 2q(/3 + M)] 

4im „ , m /o ^-2 ^ 2im , 
r*A ^ ' \ ' r 

- [2(a - /3)(r - 3A/ - 2/3) - - 3/3^ - 4Mre"^] to 



r^A 
4ir) 

-tu re 



4im 2i'-A / I o 

-e r(tj + 2t47)(l — Gs ) 



^ „ ,. „ +a) 
rA ^ ^ 

(q - 2M - f3) [2-cuC't + {Cl - 



4im^^_A 
r^A 



^i^e'-^t^' + ^e''(2tu + a;)(4Af + /3 - 3q) 
rA r^'A 



2im 



e" [2(/3 + 5a - 4,M)vj - (3/3 - 5a + 8M)c^] 



r^A 

foi7 = ^(2ix7 + c<j)(4M + 3/3 - a)C,' - ^(a + /3)ti7 



2im, 
rA 
2im 
rA 



(rA + a + 5/3 + 8M - 2r) u 



(90) 
(91) 
(92) 

(93) 

(94) 



(95) 
(96) 
(97) 
(98) 

(99) 



fo20 = ^ {4r(a + /3)ti7 + 2r^(a + 2re 



(100) 

- 2 [2a(r + ae^) + r^A] uj) Cf 



2e'' 



+ -^J [4a^ + r(l - A)a - 2r^Ae-^ - 3r^e-^ (l + 4e~^ - 2A=^)] C 



r^A 
2e^ 
r^A 
4e" 



rA H — a e + 2a 1 cij — 12(a + (5)-cli + r(a — 15re )lu 



+ — (a + /3) [2a-6r + 8Af + (A-2A*)r] n7£. 



r^A 
2e" 
r^A 



[(2a^e^ + r^A + 2ra)(A* - 1) + 2rA=^(r - 4A/ + 2a + 2/3)] ujC 
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622 



+ 



^24 
&25 
626 



627 



+ 
+ 

+ 



&28 



r- 

3e^A^ 

r5 A 

4 2 

e 

r 

46^" 



O 21^ A± 



r 2 -A / . ± 
re 



2 



+ [2(1 - 2A*)r^a; + (4e^a^ + 12ra - r^A)ro)] £ 



HA 

26" 



2(r + ae )ui — r lu \ L-^^ H -e £3 



(a + /3)(a + re"^)7i7£2 



r^A 

r^A 
_3_A^ 

r3 A 
8e^ 
rA 



[r^(A + 4 - 6A=^ + 12e"^) - 4Q^e^] a;'/: 



^"-^ ^r\-^u' + (3q - 4/3 - re-^) J A* 



(a + re 



— [2(q — 3re ^)u) — r^e "^lj'] £5*^ 

-^J^ — a — 2/3)0; — r^e^'^o;'] £^ 

e" 



+ [(6a + 3/3 - 7re"^)A=^ - 2(q + /3 - re'^)] C 



r^A 



■cj' [4a^e^ + r^(A - A± - 2 - 12e~^] £^ 



ro(a + /3) [(a - /3)e^" - 2a + /3 - re"^] £± 



(q + /3)(e'" - l)ti7 + (4 - 3A±)r-J Ci 



r^A 

e-^ 
r^A 

^e^(^-^'c.A±£± 
rA 



Background coefficients for tfie evolution equation 13211 for Hi 



Coo = 



2a^y+A 



(101) 



(A* - A) - 4rMe"^ + 2/3(r - a - 2iV/)] c^£^ (102) 



(103) 



(104) 
(105) 
(106) 



(107) 

(108) 
(109) 



coi 

C02 = 

ClO = 

Cll = 

C12 = 



= 1 



8im -A , N 

— : — e re + oi)vj 
A 

4im y 2 
— ■ — e r taj 
A 



(110) 

(111) 
(112) 
(113) 

(114) 
(115) 
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ci3 — [(2ae'^ + — r^Lj'] (116) 

ci4 = "a"'''^ ^ 
2 

C20 = A^s*^ [4(a + /3)(q + re~ )ii7 - rao;] £3 



+ (4re~^ + r(l- A^)+4a)iu'C^ 

rA ^ ' 

+ {ir^J - 2ae^a;) (118) 

C21 = — e ^ (re +0)^2 (119) 

C22 = i [(l-2A±)a;-2(a + re-^)a;']£^ + icj£j (120) 

C23 = + + + (121) 

C24 = ^e^-'^^A± (122) 

OA it 

C25 = |^e''^(a - /3)r± (123) 
Background coefficients for the evolution equation 13311 for T\ 

doi = ^[{a-Sl3){a-re-^) + C;''a{a + p)] 

+ J(2e-^-A) (124) 

do2 = -^[2re-^ +3{P-a)] (125) 

rA '- 

- ^{a + (3)e''[C;''a+{a + 2re'^)]zj (126) 

'^ii = ^e''^(/3-") (127) 

rfi2 = ^^e"^ (ae^ro — r^tj') (128) 

di3 = ^ [4(Q + /3)cc7-r(a-3re-^)a;'] (129) 

di4 = — -ime"" (130) 

dis = (131) 

die = ^ [r(a - 3re~^)tj' - 2(q + /3)ti7] - iimtj (132) 

rf20 = -jio{£t + Cf) (133) 

d2i = -^e" (2e^ati7 + rV) £^ (134) 

d22 = ^e"r'^zu£.^ (135) 

rf23 = -0(22 (136) 
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^24 = —-r-rcj'Ci' 



(137) 



Background coefficients for tfie evolution equation 13411 for hg 



eoo = 




eoi — 


6 


eio — 


i?^ r . / 
—imuj H 4( 


en = 


4im ^ 2 
— : — 6 r 
A 


ei2 = 


ij" 2 -A / 
;-r 6 UJ 

A 


ei3 = 


-612 


620 = 


4^2 /,± 

— -6 r c^7^o 
A ^ 


621 = 


A 



Background coefficients for the evolution equation 13511 for H 



koi 

ko2 

ko3 
koi 
ko5 
ko6 
ko7 
kio 
fell 



a + P 
a + P 



2r(Q + P) 



e"" (/3 - 3q - 2r-6"^) 

A6" 

(/? - 3a + 4M) 



2(a + /3) 



6~'"(3/3-a + 4M) 



a + P 



e'^iP - a + 2M) 



a + P 



S 3 h' — X 

r e 



(138) 
(139) 
(140) 

(141) 

(142) 
(143) 
(144) 
(145) 

(146) 
(147) 
(148) 
(149) 
(150) 
(151) 



a 



2im 



r^a + P) 



[{a - P){2r - 6M - P) - 2a/3 - 4rAfe"(\]52) 



— {a + P) [{Cs ~l)co + 2-coCT\ + im {cluj - O) 



rA 
2im 2 

— 6 



Ar{P + M) 
a + 



+ (a + /3-r6"^) 



+ ge-^C7? + [rA-4(a + /3)]ti7 

= ^C*? \rMl - 46"^) - 3/3 - a - 4Ml 
r^A '- 



l"l 2 

r^A 



2(4M - 3q + /3) + 



rA(8M + 3/3 - 5a) 



2(a + P) 



ki3 
ki4 



im 2 -A 
-— — C,r6 
A 



(153) 

(154) 

(155) 

(156) 
(157) 



24 Evolution equations. 



k2i = 

+ 
+ 

+ 
+ 



^ [AaCt -{a- /3)A±C?Al + 
2(a-/3) 



2f ^m('^"3M-/3)+ , 



313 2 2{P - 3q + 4M) _x 2 

~2^!i 7 — ; — S\ ^ 

r{oi + p) 



3£f + ^(a + /3)/:± 
" a + 3/3 + 2A/ ^± , A^ 



^2 -A 



^2 -A 



/oo 
/oi 

/02 
/03 
/04 
/lO 
/ll 



Background coefficients for tfie evolution equation 13611 for R 
1 



/l4 
/15 
/20 



2r4 



(/3 - 3q + 4M) 



2r2 
Q + 



-2;^(/3-3a + 4M) 



= ^-e" [(lla + 3/3-8Af)r-6Qe^(Q + /3)] 



2r4 
2im 

= -^(16M-16Q-Ar-)(x)' 



2im 

r^A 



[A(c 



r) + 8(/3 - 3M + 4q)] 



/12 = T- [rAfi + 8t^(/3 + M)] 

rA 



/13 = 



2im 
r^A 



e~''(/3-2a + 3M)tj' 



+ ^e"" [5Afr + 2ae^(2c 



r^A 
2im 
r^A 

16im 
rA 



^(/3 + 4q)] 



e"'-'(/3-M + 2Q)a; 
(M - ol)vj 



e 

r4A 

4 
r^A 



^(a + /?) [32(M - a)£± + rA/:f ] 



e"""^ (/3 - a + 3A/) cj'ZlJ 



+ 



-^e"" [2(a + /3)(rA± + 8a - 8/3) - r{M + /3)] 
-^e"''(/3 - a + 2Af) [2e^an - r'^J] Cf 



(158) 



(159) 



(160) 
(161) 
(162) 
(163) 
(164) 
(165) 

(166) 
(167) 

(168) 
(169) 
(170) 



(171) 
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/23 



+ 
+ 



/25 
/26 



+ ^e-'{a + l3){a-M)n/:t 

- ■^^e-' [a(M + f3)e^ +8{a + I3)(M - a)]ujCt 
4 A ± 

~ TsX^'" [(2°=^ + /3« - AfQ)e^ - 2r(a - M)] toC 

+ -4^^'" I" (4a - /3 + 3M)A± + 4a - 5/3 - 9Afl 
16 

4-e~""^ro [I3{a - r) - a(Q + r) + 2M(/3 + 2a)] 



2 



" [2(a - Mf + ~ I3M + 2/3a] uj'jC^ 

(a + 13) [3rMe~^ + I3{2M - r) + a(3a - 4r + 5M)] 
rj7a(a + /3)(M - a)C~^Ci 



(r — 2ae'^) i:i7£f 



|^(a-M)A±-£f 



^ (4a + /3 - 3Af ) 
^ti7(M-a)/:J 



/24 = ^(/3 + M)(/3-aV£± 



o A ± 

-^e-^(/3 + AfV£± 



4T-e"''(/3 + M) [(A* - 2)£j - 2A£f - 2£3±] - 



2 16 

— e"''(a + /3)ti7£f - ;:3^e""(a + /3)(A'/ - a)vjCt 



^ e"""^ (/3 - 2a + 3Af) aj'£^ 



(172) 

(173) 
(174) 
(175) 
(176) 



r^A 



(177) 



Background coefficients for the evolution equation l|87|l for V 

a + (3 



900 = 

goi = 
go2 = 
gio = 
+ 
+ 

gii = 



r(a + /3)e-"' 



r 

im 



-{l3 -3a + 4M) 



(178) 
(179) 
(180) 



— [2/3(r - M) - 2r(2Af - 3a)e"^ - a(2Af + 5/3) - a^] e~'''vj 



[/3(/3 - 3M) - 2(re"^ - M){l3 - a) - 4rA^e"^] 6""^ 



r^A 

im _^_A , 
— -e vj 
rA 

i^e""-'' [2(Af - a) + (a - /3 - 2Af)C5] zu 



(181) 
(182) 



Evolution equations. 



gi2 = (2^7 - Afi + ACf ro) (183) 

im _x 2 , 
qi3 = -^e r vj 

- ^ [C,2(3Q-/3 + 2re-^)+2(a-re-^)]i:a (184) 
gi4 = ^e""^ - 3q + 4Af)C,' + 2(a + 3/3 + 2M)] 

1?T? 

+ ^^e-"(3/3-a + 4M)tj (185) 



"(a^li)uj (186) 
'^Cfr^^xj (187) 



gi5 = 


2A^ 




im _ 




2A^ 




im 


gi6 = 








920 = 


A^^ 




■UJ 


921 = 













2 +^^C2 



4im^_i,_A_2 
r^A 



S20 



= -t- [(r + 8Af + 8/3 - r-C^ A=^)ro + ro;] £^ 



(188) 



(189) 



Background coefficients for the evolution equation I38II for U 

soo = ^[P{a + l3) + 2re'^(M -a)] (190) 

SOI = -^e-\a + f3) (191) 

sio = -imfi - ^ (4/3 - r + 4M) (192) 
rA 

sii = -^e-'{a + l3){4P + rA + r~2re-^)uj 



r^A 

e-"~^r\fi + M)Lo' (193) 



si2 = ^e-''-^(A.f + /3)a;' (194) 
Sl3 = — Sl2 (195) 



rA ' ' '^^^ ' -"J-^ 

+ (196) 

521 = i [r^e"^ru' + 2(Qe"^ - r)7J7 + (3a - /3 + 2re~^) ^xj] £fl97) 

522 = [(3/3-2r + 6M)a-/3' +2/3(r-3A'/) +4rAfe-^] Cj£^ 
+ [(5/3 - 6r + 14Af)a + - 2(r - M)(3 + 4rAfe-^] /I^ 

(Q + /3)£^ (198) 



rA 

-1/-A 



S23 = ^-—j^vj[{P~a + 2M)C'l+2{M~a)\C^ (199) 
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S24 = ^r'e-'-wjC^ (200) 

- ^ - 3a + 4M)C^ + (7/3 + 3a + 4M)] zuCf (201) 
S26 = — r- [(a + /3)n-(Cj + l)(3/3-Q + 4M)w]/:± (202) 



